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Abstract 

Jordan-Brans-Dicke theories with a linearized potential for the scalar field are in- 
vestigated in the framework of the stochastic approach. The fluctuations of this 
field are examined and their backreaction on the classical background is described. 
We compute the mode functions and analyze the time evolution of the variance of 
the stochastic ensemble corresponding to the full quantum scalar field in the pre- 
big-bang regime. We compute fluctuations of the term discriminating between the 
two branches of solutions present in the theory. We find, both analytically and 
upon direct integration of the stochastic equations of motion, that the dispersion of 
these fluctuations grows to achieve the magnitude of the term separating the two 
classical solutions. This means that the ensembles representing classical solutions 
which belong to different branches do overlap; this may provide a quantum mechan- 
ical realization at the level of field theory to change among solutions belonging to 
different branches. 
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1 Introduction 



In general scalar tensor theories of gravity, such as Jordan-Brans-Dicke (JBD) theories, 
the gravitational coupling constant is replaced by a dynamical field 0. The theories are 
parameterized by a dimensionless kinetic coupling parameter uj. In the limit u — > oo these 
theories go over to Einstein gravity. From time-delay experiments it is known that to must 
be larger than ~ 500 at the present epoch 0. It is conceivable that in the early universe 
the actual value of to differed significantly from that bound in being itself a function of 
the Planck mass, i. e. of the JBD field <3> = = G^ 1 (see e. g. ||). Being part of 
the dynamical evolution in the early stages $ is assumed to finally acquire a finite value 
which is consistent with observations of G and uj today. 

Theories of gravity arising from Kaluza-Klein theories or supergravity may appear as a 
JBD theory after compactification to four dimensions with uj of the order of unity with 
the square root of the metric of the compactified dimensions acting as a JBD field on the 
uncompactified space-time manifold. On the other hand, the universal part of the low 
energy effective string action || can be viewed as a JBD theory with the specific choice 
uj = —1 after a field redefinition $ = and being the dilaton of string theory. In 
this context a (nonperturbatively) generated potential for the corresponding scalar field 
is taken into account which can effectively act as a cosmological constant and drive a 
period of inflation. This has been extensively studied in the literature |J. Recently the 
possibility of having a kinetic type of inflation (in contrast to that driven by a potential 
term of the JBD field or an inflaton) has been an active field of investigation [/J. 
Especially much work has been devoted to the special case of the string low energy ef- 
fective action {u — — 1). It has been observed that there are two different cosmological 
solutions (one defined for t > to and the other for t < to), related by scale factor dual- 
ity (SFD) and separated from each other by a curvature singularity, which offer a quite 
attractive pre-big-bang scenario ||. Due to the negative kinetic coupling one of the 
pre-big-bang solutions undergoes a superinflationary period totally driven by the kinetic 
energy term. This stage is supposed to be followed after the big bang by an ordinary 
radiation dominated expanding Robertson-Walker universe. 

On tree-level, the curvature singularity prevents a smooth transition from the inflation- 
ary phase to the radiation dominated universe (graceful exit problem) ||. Later it was 
proven [ 10 1 that a dilaton potential cannot trigger such an exit. However, the singularity 
is thought to be a property of the low energy effective action which would be absent if 
higher order corrections were taken into account thus intertwining these two phases [ lT| . 
As an alternative approach quantum cosmological techniques have been applied to the 
scenario in the regime where the tree-level effective action is still considered to be a good 
approximation W*2 . 



SFD is part of a larger 0(d, <i)-symmetry |H| present in low energy action and reduces 
to SFD in the case of a homogeneous and isotropic universe. A potential for the dilaton 
as well as a curvature term break SFD. In many situations one can have an approximate 
symmetry since a spatial curvature term can safely be ignored at very early times and the 
influence of the dilaton potential is suppressed by the kinetic term as in the case of the 
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superinflationary branch close to the curvature singularity. General scalar tensor theories 
such as the JBD theory don't possess this kind of 0(d, <i)-symmetry which is clearly of 
stringy origin. Nevertheless solutions can be found which form two disconnected branches 
separated by a field dependent discriminant; they are (apart from singularities) defined 
for all times such that one always finds similar pre-big-bang scenarios as those described 
above for the stringy case, but one also faces similar problems. In JBD theories, i. e. 
those with arbitrary but constant u, such periods of kinetic inflation are present only if 
the parameter u is negative 0. These models will be the most interesting ones for our 
purpose 
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Quantum fluctuations in JBD theory have been studied on inflationary backgrounds, 
i. e. those fulfilling the slow-roll conditions, in | I5[] , where also density perturbations 



have been computed; similar analyses have been performed on more general Robertson- 
Walker backgrounds, but without potential, or alternatively for an exponential potential 
in the context of extended inflation in the Einstein frame [16) . For the string cosmological 
scenario (cu = —1) scalar and tensor perturbations have been computed |T7| and the exit 



problem via quantum backreaction has been investigated in |18| . We point out that a full 
analysis of the backreaction problem in JBD theory on general cosmological backgrounds 
including a generic potential (and also a matter contribution to the energy density) is still 
lacking. One step into that direction shall be the subject of this paper. For the problem 
to be tractable analytically one could e. g. expand the solutions for small values of the 
potential or for small times where the influence of the potential can be linearized. 

The intention of this paper is to study the dynamics of quantum fluctuations of the JBD 
field and its backreaction on the background fields by means of a stochastic analysis. For 
this purpose we split the field and its first time derivative (which will be treated as an in- 
dependent stochastic variable) into sub- and super-horizon sized modes. This technique 
has been applied to similar problems before fH| . In this framework it is possible to trace 
the time evolution of the long wavelength (coarse grained) fields which we interpret as the 
classical background under the influence of the stochastic noise operators. It is possible 
to write down a system of stochastic differential equations (Langevin equations) and the 
corresponding Fokker-Planck equation for the probability distribution of the stochastic 
variables, i. e. the classical background field and its first time derivative. The Langevin 
equations can be integrated numerically by a well known method [2C, For the Fokker- 



Planck equation one can in principle search for an analytic solution. 
The stochastic method is suitable to investigate the effects of quantum effects on the 
background quantities as long as the size of the fluctuations are under control. It is a well 
defined technique and can serve as an alternative to the minisuperspace approach of ex- 
amining the exit problem of string cosmologies. Since the stochastic differential equations 
can be integrated for any form of the potential one can hope to make statements about the 
influence of the scalar field potential in this context. However, we argue that a linearized 
potential (which shows up as a cosmological term in the action) will be sufficient as a 
first approximation. For the linear potential we will also investigate the mode functions 
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on that background and their dependence on the potential and the parameter 00 since 
they will finally enter the stochastic differential equations. By comparing the relative size 
of the mode functions and the zero mode it is also possible to check the validity of the 
procedure. 

With the above tools we are able to study the backreaction of the fluctuations on the 
dynamical level. We check up on the role of a large potential and examine how it decou- 
ples from the dynamics in the limit when it becomes small. We compute the variance of 
the scalar field fluctuations and its backreaction on the Hubble parameter which we will 
finally compare with its discriminant separating the two different background solutions. 
If both are approximately equal, it is no longer possible to unambiguously talk about two 
distinct branches. 



In section |] we give the basic equations of JBD cosmology and fix our conventions. The 
setup of the stochastic treatment of the scalar field fluctuations are presented in section ||. 
We discuss the general background solutions for a linear potential in JBD cosmology and 
give expressions for the mode functions on that background in section |] together with 
numerical solutions for variable uj and strengths of the potential; we analyze the exit 
problem of the superinflationary branch in the presence of quantum fluctuations of a 
scalar field. The results are summarized in section 



2 Basics of Jordan— Brans— Dicke cosmology 

We start with the following action of the JBD theory, 
1 
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including matter and a finite potential energy term for the JBD scalar field $. TZ is the 
Ricci scalar of the metric. In the following we keep u as an arbitrary parameter of the 
theory assuming that it is independent of the field $. Since the physics is more obvious 
in the JBD frame, we will not change into the conformally related Einstein frame. 
Varying the action ([[]) with respect to and $ yields the following two equations of 
motion 
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where T mfll/ is the energy-momentum tensor of the matter sector. Combining the trace 
of equation (|2]) with equation (|3]) one arrives at a Klein-Gordon equation for the JBD 
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scalar $ 



n* = -*!-T m —*-v W+ 1 (4) 

2u + 3 2uo + 3 V ; 2uo + 3 3$ ' V ; 

with T m = being the trace of the energy-momentum tensor. For the homogeneous 
and isotropic metric we use an ansatz of the Robertson- Walker type defined by 

g^ dx^dx u = dt 2 - a 2 {t) dQ 2 3 {k) . (5) 

Here a(t) is the scale factor of the universe and dQ 2 (k) stands for the volume element of 
the 3-dimensional spacelike hypersurface of constant t and k — ±1,0. In this metric an 
ideal fluid is described in its rest frame by T m = p — 3p (with p and p denoting the matter 
energy density and pressure respectively). The equations (§) and ((^) in the metric @ 
consist of only two linearly independent equations, the third one serving as a constraint. 
We choose the 0-0 component of (|2|) in the metric (|^) which reads 

""3$ 6 $ 2 $ 6$ 6 a 2 $ 2 3 a 2 $ a 2 ' 1 ' 

where we have used the definition of the Hubble parameter, H = a/a and an overdot 
denoting a time derivative. For later purpose it will be useful to write (|6]) as an explicit 
equation for the Hubble parameter, 




2$ V 3$ 12 $ 2 6$ 6a 2 $ 2 3a 2 $ a 2 ' K} 

which shows that there are two different background solutions separated by a discriminant 
given by the square root in (|7|). By setting uj = —1 for a homogeneous and free field in 
a universe described by p = k = it becomes apparent that the + sign (— sign) in (|7|) 
corresponds to the (— ) branch ((+) branch) of reference |§ defined only for t > (t < 0). 
Therefore, the discriminant distinguishes between a superinflationary and a collapsing 
phase of the (+) branch for negative times and between an ordinary radiation dominated 
expanding universe and a collapsing one of the (— ) branch for positive times. At each 
instant of time the separation originating from the square root of (0) is, however, equal 
to the separation between the two different branches if they were defined for all times. In 
subsection |4.3| we shall compare the discriminant with the fluctuations in the JBD field. 
Equivalently, the Klein-Gordon equation @ in the Robert son- Walker metric (|5|) is given 

by 

a 2 2w + 3 2w + 3 1 ; 2w + 3 W K J 

The system of equations (fj)-© contains also spatial inhomogeneities of the scalar field 
$(x, t) which we assume to be small but non negligible and living in a homogeneous 
universe described by the scale factor a(t). Since the purpose of this paper is to study 
the behaviour of inhomogeneous modes of $ we separate them from the homogeneous 
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background, i. e. the zero mode of $. We decompose the scalar field $ into a homogeneous 
background and an inhomogeneous part according to 



$(x,t) = $ (t)+5(x,t) , (9) 
and expand the latter within a box of volume v = I 3 into its modes as follows : 

<5(x, t) = T 3/2 £ [a k <p k (i) e* kx + a], ^ (t) e^} . (10) 

k 

The sum is performed over all wave vectors k = ^ with n 6 Z 3 . The equation 
of motion for the homogeneous zero mode $o(£) is clearly given by (§) with the spatial 
derivative dropped. On the other hand, substituting the ansatz (H) into the Klein-Gordon 
equation expanding to linear order in <5(x, £) and using the equations of motion for the 
background $o(X) we immediately arrive at the equation of motion for the mode functions 



y5 k + 3Hip k + 
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2u + 3 



V'(%) - V"(%) $ ( 
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with V'($o) denoting 9V^($)/9$|$=$ , etc. This equation is independent of the curva- 
ture term k/a 2 which means that the topology of the spatial section enters only via the 
background quantities. 

For the rest of the paper we will neglect the spatial curvature term (k = 0) as well as 
the contribution of ordinary matter to the energy density (p = p = 0). For completeness 
we give the background solutions in the case where in addition the potential is neglected. 
The equations of motion @ and @ allow the following power-law solutions which are 
defined for all times except t = 0:Q 

t\ riuj) and aft) ~ Itl'M , 



$o (t) 



with the exponents given by 
r{uj) 



1 ± + 



4 + 3cj 

After the field redefinition $ 
action (uj = — 1) |§ 



and 



(12) 



(13) 



e * we recover the known result for the low energy string 



-rJ_ 

(f)(t) ~ (^^3 - 1) In 1 1 1 and a(t) ( 14 ) 
in which case it is invariant under SFD. Comparing ( flljf) with (|14|) it is clear that the 



+ sign in front of the square root in (13) corresponds to the accelerated expansion in 
pre-big-bang cosmology, i. e. the superinflationary part of the (+) branch. For t > the 
part of the (— ) branch which corresponds to the — sign in (|T3|) can be smoothly connected 



1 the singularity can certainly also be at t = t - For simplicity we set to = 0. 
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to a FRW universe with decelerated expansion and can therefore account for our present 
universe. We note that there is also a constant background solution for the vanishing 
potential, i. e. $o = H — 0. We will not be interested in that possibility, however, since 
it implies a static universe. 

For the zero potential background the mode functions can be directly read off from equa- 
tion (|TT|) and given in terms of Hankel functions, 

^ k (t)oc|t|-^c H (/3Mitr^) , (15) 

where for shorthand we introduced the quantities 

/ s r(uj) „. , 1 Ikl , , , . . . 1 1 — 3q(u) , , 

a( U ) = ~U, m = T -^A 7M = 1- 5H , Ku ) = ± -_«U. (16) 

By ao we denote any initial value for the scale factor which we can set to one for simplicity. 
C u stands for any linear combination of the Hankel functions and with the two 
integration constants Ci(k) and C2(k) multiplying the respective Hankel functions. Choos- 
ing a certain linear combination corresponds to making a specific choice for the vacuum 
state; it is subject to the Klein-Gordon normalization condition |c2 (k) | 2 — |ci(k)| 2 = 1. 
The correct Minkowski limit (positive frequencies) can be obtained by the simplest choice 
c\ = and C2 = 1 for all wavenumbers which is the vacuum state we shall adopt in the 
following. Choosing the negative sign in the above expression for gives the non triv- 
ial solution $ ~ 1 1 \ r ^ in the limit k — > 0, whereas the positive sign yields the constant 
solution. 



3 Stochastic equations and correlation functions 

In this section we will derive the modified equations of motion in the stochastic approach 
and the correlation functions between the various stochastic field operators. For this 
purpose we split the field and its first time derivative in sub- and super-horizon sized 
parts |3 [HJ, 

$(x,t) = $<(x,t) + $>(x,t) 

(17) 

$(x,t) = v<(x,t) + u>(x,t) , 
and decompose each of them into Fourier modes according to 

$ > (x, t) = r 3/2 J2 0[± ( k - eaH )\ («k Mt) e Jk ' x + a{ iflit) e" lk - x ) (18) 

k 

v> (x, t) = r 3/2 J2 ( k - ™ H )\ («k Mt) e Jk ' x + 4 ^k(*) e-**) • (19) 

k 

The quantities $< and v< are the long- wavelength (coarse-grained) fields and correspond 
to the lower sign in the 9 functions of equation fll8|) and (0), whereas $> and v> are the 
short-wavelength fluctuations (referring to the upper sign). The parameter e denotes the 
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splitting point and should be of the order of unity if we choose the horizon length as our 
effective physical cutoff. We define the stochastic noise operators 

X(x,t) = ea(H 2 + H) r 3/2 ^S(k-eaH) (a k <^ k (t) e ik ' x + aj, ^(t) e" ik - x ) 

k 

£(x,t) = ea(H 2 + H) r 3/2 J2 S ( k ~ £aH ) (a k ^k(t) e ik ' x + aj, k (t) e^ k ' x ) . (20) 

k 

In the following we will aim at deriving a set of stochastic differential equations (Langevin 
equations) describing the system given by the action (H). We regard the short-wavelength 
quantities $>(x, t) and f>(x, t) as small perturbations around the coarse-grained fields 
$< (x, t) and t>< (x, t) which we will finally interpret as semiclassical random variables 
whose evolution we follow under the influence of the stochastic noise terms (|20|). It is 
easy to see that 

$<(x,t) =u<(x,i)+xCM) , (21) 
where we have used the relation 

Q 

— 9[ =p(/c - saH)\ = ± ea(H 2 + H) 8(k - eaH) . (22) 

In the same way we get $>(x, t) = v > {x, t)—\ (x, t) and hence $(x, t) = f>(x, t)+f<(x, t). 
Equation ( pi]) is the first stochastic differential equation since it contains only coarse- 
grained fields (our random variables) and the stochastic noise. 

It is more difficult to obtain a Langevin equation for f<(x, t). We write ti<(x, t) = 
$(x, t) — v > (x, t) and get 



t)<(x,f) = -3ff6(x, t) + A ^ t} - ( 2y(d>) - $(x, t) 



(23) 



+ £ (x, i) - V 3 ' 2 °( k - £aH ) k e ik ' x + 4 <p* k {t) e 



ikx 



where we used the equation of motion for the full (inhomogeneous) quantum field (||) and 
the definition of f>(x, t), equation (|T9"D. The first line of equation (|23| ) can be split into 
the coarse-grained and short wavelength fields according to (|T7|) . We will expand the 
potential and its derivative around the coarse-grained fields, 

2V($) - V'(<$>) $ = 2V - V ' $< + [Vq- Vq $<]$> + C($>) , (24) 

where we have used the notation V = V(§ ), V ' = dV($)/d&\$ = $ , etc. To simplify the 
expressions we will approximate $o ~ and vq ~ v < which means that effectively we 
neglect all spatial gradients of the coarse-grained field $<. Later we will see that this is a 
very good approximation; it gets exact in the limit e — > 0. This identification will cancel 
the terms linear in the small quantities $> and v > in equation ([23|). In the second line 
of (E3f) the second time derivatives of the mode functions are eliminated via their equation 
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of motion (IT). 

This together with the approximation mentioned above gives the following differential 
equation for the coarse-grained field t><, 



(x,t) = -3Hv<(x,t) 



2u + 3 



2V -V^ < (x,t) + £(x,t) 



(25) 



up to second order in the small quantities. 

The equations (|2TJ) and Q2"5l) constitute a system of stochastic differential equations for 
the fields $<(x, t) and f<(x, t) which we regard as independent random variables subject 
to the two independent stochastic forces x( x , t) and £(x, t). Effects of spatial variability 
of the fields enter only through the stochastic noise operators. For a full description of 
the system an equation similar to (0) should be added to the two Langevin equations; the 
total energy density we assume to be stored mainly in the background fields. Therefore, 
we neglect energy contributions arising from the small quantities $>(x, t) and f>(x, t), 



2 $^ 



2u + 3 v% 



\ 12 



6$< 



(26) 



The set of equations (fZTJ), (j2"5| ) and fl2fi|) can be integrated by modeling the stochastic 
noise terms by a random number generator. Statistical independence of the numbers 
requires a Gaussian distribution. Integrating Langevin equations of that type is a well 
known problem and one can use standard techniques . 

The correlation functions between the stochastic noise operators can be calculated in a 
straightforward way, 



(0| X (x,t)x(x',t / )|0) 



1 £ wfl + £l sin( g qg|x-x'|) x 



2tt 



\k\=ea\H\ 




eaHlx — x' 



(27) 



and similar expressions for the correlation functions ( | £ (x, t) £ (x', t') | ) and the cross 
correlation functions ( | x ( x , t) £ (x', if) | ) and ( | £ (x, t) x (x', t') | ), for which |v?k(^)| 2 
is replaced by |0k(£)| 2 , Vk(i) ^if) an d y?k(£) y?kW respectively In the same way one de- 
rives the only non trivial commutator for the noise operators, 



[x(x,t), £(x',t' 



x 



(2* 



|k|=eo|H| 



8(t - 1') . 
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4 Mode functions with a cosmological term 



The aim of the following section is to analyze the mode functions of the JBD field on a 
background evolving under a linearized version of the potential V(<&) in the action (|I|). 
This linear term can be interpreted as a non vanishing cosmological constant, or as a mass 
term for the canonical mass dimension 1 field \ defined through x 2 = $. We choose a 
linear potential first of all because it is tractable analytically. Moreover, since our main 
point of interest is the superinflationary branch before the singularity, it is not the minima 
of the potential, which are modeled by nonlinear terms, which control inflation. Actually, 
the influence of the potential on the dynamics is naively expected to be subdominant, 
because, for small values of u, the grander part of the energy is stored in the kinetic term 
of the scalar field. However, the strength of the potential can affect the rate of evolution 
which in turn can change the magnitude of quantum (represented here as stochastic) 
fluctuations. To study this effect in the system not too far away from the singularity the 
linearized version of the potential term is expected to serve sufficiently well. The mode 
functions enter e. g. in the correlation functions (p7|) and the equations for the variance of 
the fluctuations. Background solutions for a finite cosmological constant have been given 
for the first time in [^2] and applied to pre-big-bang cosmologies in [ 14 ] . For completeness 
we give a short review of the derivation in the following subsection. 

4.1 Background solutions for non zero cosmological constant 

Using the ansatz V(Q) = — 2A$ for the linearized potential the 0-0 component and 
the i-i component of the Einstein equation (^j and the Klein-Gordon equation (|j) read 
respectively (spatial curvature and matter contributions are omitted) 

3F 2 -A = |^ 2 -3#^ (29) 

- 2 *-H 2 + \ = ^ 2 + i°+2i7tf (30) 
a 2 $ 

— + = — ^— , (31) 

$o 2cj + 3 ' V ; 

where we have used the definitions H = a/a and \I/ = $o/^o- Defining the following 
functions 

f = %a 3 h = %a 3 g = %a 3 (32) 

and noting that a/a = H + if 2 and $o/^o = ^ + we nn d that the linear combination 
of equations ( p9[) - (|30|) - ~ (|31|) reduces to the familiar form 



/ - b 2 f = where b 2 = 2A ^±1 . (33) 

2uj + 3 
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The general solution of this equation is / = A + e +A- e 
A + and A_ to be determined later. The lhs of equation 
can be written as 



~ bt with constants of integration 
ill) is equal to h/f so that it 



h-d 2 f = 



where 



d 2 



2A 



2u + 3 



(34) 



with the first integral h = C + d 2 J fdt (C being an arbitrary constant). From fl32|) we see 
that f — h + 3Hf which can be solved for H and since the functions / and h are already 
known it can be integrated to yield the scale factor a, 



H 



f-h 
3/ 



and a = a \ f | 3 exp 



dt 



(35) 



From ( p2[ ) we can read off that \I/ = $0/^0 — h/f and by integration we get 



$0 = $o 0> ex P 



dt 



(36) 



In general the Hubble parameter H will have a singularity. Without loss of generality we 
can fix it to be at t = 0; this requires to choose A + = —A- = ±A with A > 0. Then 
/ = ±2Asmh(bt) and the sign ambiguity determines the two branches of solutions. It 
turns out that this fixes the ratio C/A as well meaning that the solutions for H, a and 
$0 that we get by successive integration of the above equations are only fulfilled if C/A 
takes the value 



C 
A 



6A 



3w + 4 



We fix the constants of integration do and $0°' by making contact with the well 
results for u = — 1 in the limit A — > 0, a = a \ t | ±:l /v3 anc [ Q = — (1 =p In 
being again the stringy dilaton, O = — ln$ - 
The relevant background quantities parameterized by uo and A are 



H(t) 



a(t) 



$o(t) 



2A(w + l) 



V2cj + 3 y/3u + 4 



coth(fct) =F 



a 2 T 9^+12 I ft 



sinh(6t) 




sinh 1 (bt) 



9w+12 



1 



1 v / aJ+9 



sinh(fet) 



3a;+4 




. V6cj+9 



3a;+4 



(37) 
known 

1*1, 00 



(3? 



(39) 



(40) 
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where we have used q{oS) and r(oS) as defined in (0). Choosing uj = — 1 these equations 
reduce to 



Hit) 
a{t) 



=F 



^ sinh- 1 (V2At) 



1 1 

* 



in J_ 



$ (0) 2 ±v / 3 



±4 



tanb [ - V2Xt 



a 



-1=FV3 



sinh(V2At) 



1 / \ ±V3 

tank ( - V2At 



(41) 
(42) 
(43) 



where the expressions after the arrow are valid in the limit y/\t <C 1. The prefactors in 
the equations ( |39D and ( |40| ) are essentially the constants a and $q 0) and have been chosen 
such that in the special case w = -1, A = we get the simplest possible initial values. All 
the upper signs correspond collectively to the (+) branch showing accelerated expansion 
for t < and the lower signs to the (— ) branch with decelerated expansion for t > 0. 
One comment is in order concerning the sign of the potential. Up to now it was assumed 
that A is positive; this is indeed not necessary: If A is negative the correct expressions are 
obtained by substituting £» — ^ | £> | or A — ^ | A | respectively and simultaneously replacing all 
hyperbolic functions by their trigonometric counterparts. 



4.2 Mode functions for non zero cosmological constant 

Being equipped with the full expressions for the background with linearized potential we 
are now in a position to investigate the influence of this potential on the mode functions 



determined by equation (|11|) . It possible to examine the dependence on the JBD param- 
eter u) equally well. It is clear that the equation of motion for the mode functions in the 
background fl38|) - (|1]) cannot be integrated analytically in full generality. However, it 
is possible to examine the system in the limit \f\t 1. This can either be the small 
potential limit for unconstrained time parameter or an expansion around t = valid for 
general A. In the small A limit it can be studied how the potential effectively decouples 
from the solutions. Finally we will perform a numerical integration of the most general 
equation (|TT|) and compare the results with the approximate formulae. We will later in- 
vestigate the region of validity of the stochastic approach. 



First we want to give the asymptotic expressions of the solutions to (pTf ) in the limit 
bt ^> 1 for positive A (this is equivalent to \f\t ^> 1 for all values of uj we are interested 
in); in this limit the functions H(t) and a{t) in (|38|) and (|3^) can be simplified significantly. 
We find one exponentially growing and one decaying mode, ip^ ~ e Ult and ip^ ~ e n2t 
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with the constants Q\ and £l 2 given by 



3V2A(a; + l 



± 



2 I V2o; + 3 V3w + 4 \ 



5A + 



18A (cu + 1) 2 
(2a; + 3)(3a; + 4) 



The two constants Q\ and ^2 are the same for &o£/i branches. For u 
f2i ; 2 = ±v2A; for that case we also find a correction, 



± 



\ 



2A 



«o 



2 X± ^A lT ^3 



(44) 

•1 this reduces to 
(45) 



where the ± under the root distinguishes the (+) branch (upper signs) from the (— ) 
branch (lower signs) and the ± in front of the root refers to the two different constants 
Qi and ^2- We note that similar asymptotic solutions can not be given for negative A 
because the trigonometric functions are undetermined at infinity. 

Next we can try to solve equation QTTT) for small arguments b t <C 1 where the hyperbolic 
and trigonometric functions show the same behaviour. The approximate equation of 
motion for the mode functions now read 



3?M . 

y?k h — ryr^ k + 



t \-2q(u>) _ 



2A 



2c; + 3 



y?k = 



(46) 



and is valid for all uj not too close to the singular point —3/2; an overdot denotes a 
derivative with respect to \ t\ . The dependence on A and its sign enters only through the 
last term in the brackets. It is clear that the solutions only depend on the modulus of t. 
For k = we obtain an expression which passes over to the zero mode in that limit and 
for A = we recover the free solution (|15|). 

Unfortunately this equation is still too difficult to solve in general the trouble coming 
from the irrational powers of the first term in the bracket; it cannot even be solved as a 
power law series. However, for the most interesting case of uj = — 1 on the (+) branch 
the exponent can well be approximated by — 2q^ (uj = —1) = 2/y/3 ~ 1.15. So one can 
look for a solution with the exponents 1 and 2 and compare the two cases. Later we will 
comment on the validity of this simplification. We want to solve for uj = —1 the equation 



V3. 



t 



2A 



= , 

with m replacing 1 or 2. This equation can be solved via the power law series 

oo 



(47) 



(4* 



n=0 
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The exponent r and the coefficients c n can be computed in the standard way. For m = 1 
we get the two linearly independent solutions 



c 



1-±±^AI«I 



\/3|k| 



a 2 , 



t\ 



2V3 



12 



A 2 |t| 4 - 



23 + 9^ |k| 2 



390 



a 2 



Alt! 5 --- 



(49) 



c 1 1 



1 1+^3 



v/3 



Alt 



4- v^lkl 2 



a 2 , 



39 

381 - 157^ Ikl 2 



132 



12870 



a 2 , 



AU| 5 + - 



(50) 



and similarly for m = 2 

*£>(*> 



Co 



l-i±^A ltl 



3 + 2^ 



2A 2 



21 + 13^3 



1584 



24 

2A 2 - (7-3^) 



1 



- V3) 
Altl 6 - 



(51) 



c 1 1 



a+V3 



, 3- \/3 , , l2 9-4^ 
1 + — —A t 2 + 



+ 



6 

75 - 37\/3 
36432 



264 

2A 2 - (13 + 3^) 



2A 2 - (3 + V3) 



Altl 6 -- 



(52) 



These are two linearly independent solutions for each case with two different arbitrary 
constants Co and Co parameterising them. Differences between the solutions for m — 1 
and m = 2 are of third order in 1 1 | . It can be read off that the zero mode in the free case 
(A = 0) is a constant or goes to zero as ~ 1 1 | 1+v/ ^ in accordance with what we found in the 
previous discussion. Comparing with equation ( fj3|) the two yet undetermined constants 
are fixed to be c = c = $ 0) . We are not interested in the constant solution because we 
want to avoid a static universe. 

Recursively, the higher order coefficients can be computed according to the following 
relations: 

11,12 



2Ac 



n-2 



C n -3 



^2 j /y 



(2n - 1 - \/3) + n (n — 1 — y/3) 



for m 



(53) 



2Ac 



n-2 



Cn— 4 



r 2 + r (2n - 1 - \/3) + n (n - 1 - VS) 



for m = 2 . 



In figure |T| the results of a numerical integration of equation (|i~T|) are shown for = — 1 
and the linear potential discussed in this section. The qualitative picture is the same, 
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Figure 1: Mode functions for u — — 1. In (a) they are shown for various wavenumbers 
|k| as indicated in the figure for negative times. The potential was set to A = 1. The 
zero mode corresponds to the bold line. In (b) the exact solution is compared with the 
approximate ones given in equation (^) and (|52D for A = 1 and |k| = 10. Accordance is 
achieved for small t. 

however, for general u: Starting from a power law growth for the zero mode as \ t\ is 
getting larger, this behaviour is converted into oscillations with increasing amplitude and 
frequency after a characteristic time for the mode functions with larger k. This initial 
growth is well approximated by the power law series solutions as can be seen in (b) where 
the full untruncated solutions (|50|) and ( |52l) are shown. The mode functions are similar for 
all k around the origin and ^k/^o — * 1 as one approaches zero. We point out that for the 
linear potential it is not necessary that the higher mode functions are much smaller than 
the zero mode, because in this case equation ( |XT| ) is exact and there are no corrections 
cx |y?k| 2 etc. 

Away from u = — 1 we shall be contented with the qualitative picture given by the 
numerical integration shown in figure |2|. Although u is a free parameter of the theory 
we are especially interested in u lying around —1. Higher dimensional theories in D 
dimensions compactified to a maximally symmetric internal manifold yielding a d + 1 
dimensional Kaluza-Klein cosmology can be viewed as a JBD theory described by a 

parameter u = — 1 + j=>; the role of the JBD field is played by \J g ( mn [|]. Going from 
uj = —1 towards somewhat larger values accessible to the above formula it is apparent 
that the amplitude of the mode functions are significantly suppressed. They become 
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Figure 2: Mode functions for u « — 1. Their dependence on is demonstrated. A 
sample of solutions for tp^ is depicted in (a) where A = 1 and |k| = 20. The values of uo 
vary between —1.05 and —0.95 the bold line denoting uo = —1. In (b) A = 1 and |k| = 1 
has been chosen. 

enormously large if u> drops sightly below —1. 
4.3 Scalar field fluctuations 

The presence of two different solutions, a generic feature of all JBD theories, have their 
origin in the sign ambiguity in front of the square root of equation ([/]). As long as only 
the classical theory is considered the root is a well defined discriminant separating the 
two branches. Fluctuations in the JBD field can, via backreactions, influence the ho- 
mogeneous degrees of freedom such that their behaviour has to be described in terms of 
stochastic quantities as well (e. g. mean value, variance). As soon as the induced fluc- 
tuations in H are comparable in size with the root distinguishing the two classical paths 
it is conceivable that the system undergoes a transition from one trajectory to the other 
and a strict distinction between the two trajectories present in the classical picture is no 
longer justified. 

The subject of this section is to make this idea more precise. Let us begin with a vanishing 
potential. According to our choice of the vacuum the mode functions are given by (p^ = 
\t\~ a H^(P\ 1 1 7 ) together with equation (^) in which we shall always choose the negative 
sign for v in order to deal with the cosmologically interesting solutions only; in this case 
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a + 7z/ = 



— r. The following relations will turn out to be useful: 
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As was demonstrated e. g. in the variance in the fields can be given, due to the 
stochastic differential equations, by the correlation functions of the corresponding noise 
operators. They are given by equation fl27D where we shall adopt the specific choice e — 1 
2 . For simplicity we neglect the spatial dependence of the correlation functions (i. e. we 
consider local fluctuations only) and define the variance of the stochastic field operators 
to be the following positive definite expressions 



(5$ 2 ) 
(8v 2 ) 
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(55) 



where the cross terms have been defined in this way in order to avoid complex quantities. 
Using equation (]27D together with the given mode functions we finally get with the help 
of (M) 
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2r 
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1-9 
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(56) 



2 In an earlier analysis the semiclassical limit e — > had to be chosen because it is only in this 
limit that the mode functions can be approximated by those of a massless free scalar field. Since we don't 
use such an approximation here we need not perform the limit e — > 0. 
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(6$6v) = (5<5> 2 )\H 
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For the last equation the commutator (^) was used. The explicit time dependence is 
reexpressed in terms of the background variables, equation (|12|). The above equations 
depend on an overall parameter a® which we shall take to be of the order of unity. 
An especially interesting example is the pre-big-bang scenario of |J which is modeled by 
the choice ui = —1. In the following we will consider the superinflationary (+) branch 
realized by q = —l/y/3 and negative t. As can be seen from the first line of the above 
system together with the time dependence of the background quantities, the variance of 
$ grows for the (+) branch like (— 1)~ 3 as one approaches the curvature singularity at 
t = 0. The additional factors of H and H 2 in the expressions for (5<&5v) and (5v 2 ) 
respectively make those variances grow even faster, i. e. (5<&5v) oc (— 1)~ 4 and (Sv 2 ) oc 
(— 1)~ 5 because the Hubble parameter scales like —q/(—t). Hence, the growth of the 
fluctuations is not bounded from above and one can easily imagine them dominating over 
the classical evolution. In particular, if 5H rms = ( 5H 2 ) 1 / 2 , the fluctuation of the Hubble 
parameter induced via backreactions, is of the same order of magnitude as the square root 
of equation ([!]) there is no need to assume that the field will stay on its classical path. At 
this point the field can be imagined to be part of either branch. 

Next we will derive a quantitative estimate for the occurrence of this effect. The quantity 
which distinguishes between the two branches is the expression under the square root of 
equation fl7|) which we will denote as A, 
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2u + 3 6§ y($ c 
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<J> 2 



6$r 



(57) 



Note that A is completely given by the homogeneous field; we are looking for its variation 
induced by substituting $ — > $o + 8$ and v o — > vq + Sv. The variation of A to second 
order gives 
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(58) 
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Since we neglected the spatial dependence of the correlation functions, the above variation 
is due to local fluctuations of $ and v and it is consistent to consider the expression (0) 
with the spatial derivatives dropped. The variance of the quantity A is 
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(59) 



We shall identify (<5$ 2 ), (Sv 2 ) and (<5<Mt>) with the expressions given in equation fl56|) . 
As a condition for an overlap of the two distinct classical trajectories induced by stochastic 
effects we require 5A rms /A ~ 1. For general u and a vanishing potential we get 



5A r 



A 



\ 



(<5$ 2 ) t (5v 2 ) n (8$5v) 

$o^o 
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(60) 



In the model considered so far 5A rms /A ~ 5H rms /H (apart from a trivial factor 2 coming 
from the square root) which means that right at the point when an overlap becomes 
significant a large fraction of the total energy is stored in fluctuations and spatial gradients 
may no longer be neglected. Although no assumption is necessary concerning the ratio 
ific/ $o for a vanishing or a linearized version of the potential the energy density due to the 
fluctuations should not be dominating over the contribution coming from the background. 
From this point of view the condition 5A rms /A ~ 1 will be met at the border of validity 
of this method. 

Setting w = -1 we find a quantitative estimate by inserting the expressions derived in 
equation (|56[) into (|6( 

S A 
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37r J3 + V3 $o 



(61) 



7.49. 



with the index of the Hankel function given by u = ~f f^7§ an d ( t^t^ 
In the second line we have used the time dependence of the background fields for the pre- 
big-bang branch. The result depends on the initial values ao and which we assume 
to be of the order of one. If we had not made the choice e = 1 at the beginning the result 
would also depend on e 3 ^ 2 which can always be absorbed by the initial values. No matter 
what the prefactors are, the strong time dependence with the power of-§ - -4.23 
ensures that the condition SA rms /A ~ 1 is unavoidably met at a time which can be read off 
from equation (|6"lD ; this time is ~ 1 if the initial condition cCq 2 / is of the order of one. 



The same as was demonstrated above for the stringy example oo = —1 can very eas- 
ily be generalized for arbitrary u (again restricting ourselves to a vanishing potential). 
Equation ( |B"UD is valid for general u. The structure of 5A rms /A is very similar to that 
of®, 

F(u) (tfai*^) , (62) 
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with F(u) denoting a function of u which is approximately one in the region of interest. 
Using the general background solutions for zero potential ( |12D and (0) it is obvious that 
5A rms /A ~ (— t)~ 2+ ^ q ^'~ 2 r H for the (+) branch. The exponent, being a function of u, 
is always negative (as it should be) and a monotonically growing function increasing from 
— | — v^3 at u — — 1 to — | at cu = 0; this is again a sufficiently fast growing function 
for the value of SA rms /A becoming one at a time which is of the order of unity. The u - 
dependent prefactor of equation ( |5^) shows the remarkable feature that it has a singular 
point at uj = where it becomes plus infinity. This means that approaching u = from 
either direction drastically enhances the ratio SA rms /A such that the overlap condition is 
easily fulfilled. 

In trying to explore the influence of a linearized potential on the ratio 5A rms /A one 
would have to refer to the approximate mode functions we found in (|50| ) and (|52|) for 
<jj — — 1 and insert them into the corresponding equation which is similar to flBUl) . ^From 
the resulting formulae it is, however, hard to see how a small but non zero potential 
influences the effect described above. Therefore, we shall rather directly integrate the 
coupled stochastic equations of motion (^1|), fl25|) and ( |26| ) using the approximate mode 



functions valid for a/A£ < 1. The discretized version of the equations can be written in 
the form 



I t n "T V6i w<p (63) 
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v n+ i = v n -3H n v n 5t + 2\<S> n 6t+^- ( 1+ * 3)2 |Cl«o 2 * 



1 + (^4^vl)^'»' 2 )'^r^^^ (64) 



TT + 
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with n being the step number and <5t the step size. To simulate a discrete Wiener process 
we use the two independent random Gaussian deviates and of unit variance and 
zero mean, 

= and (^^> = ^ mn . (66) 



The above discretized equations are in agreement with Ito's definition |23| of the stochas- 
tic noise. It is assumed that other interpretations of the noise term have only small effects 
on the qualitative behaviour of the system The Gaussian deviates are produced by 



a standard method |24[ from a very well tested random number generator. 
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As an example, the results of the numerical integration of the equations (^3) - fl65] ) are 
shown in figure |3|. The integration was done for the (+) branch (referring to the lower sign 
in equation (pop) and negative t. Both plots show two classical trajectories (obtained by 
switching off the stochastic noise terms) corresponding to accelerated expansion (H > 0) 
and accelerated contraction (H < 0). On the classical level there is no transition possible 
between the regions of positive and negative H. In both examples the initial conditions 
of all stochastic trajectories are chosen to emerge from the same starting point which is 
also the starting point of the classical path for H > 0. In (a) the free parameter |$[) 0) | a^ 2 
is varied which effectively sets the strength of the fluctuations. While smoothly turning 
on the fluctuations by increasing the value of |$o 0) | a o^ 2 the stochastic trajectory bends 
towards the second classical path characterized by negative H. Phenomenologically this 
situation is of course not attractive; on the other hand this provides an example that via 
backreactions two a priori disconnected regions in the space of solutions to the classical 
equations of motion are no longer separated when sizable fluctuations enter the analysis. 
Certainly not all stochastic trajectories behave as the one depicted in (a). In plot (b) 
the path labelled by "5.0" in (a) is redrawn together with other realizations, i. e. all 
the parameters are held fixed but the initialization of the random number generator is 
different. The dispersion of the whole ensemble of stochastic trajectories grows much faster 
than the difference between the two distinct classical solutions; this is in accordance with 
what we found analytically in this subsection. 

The condition \f\t <d 1 is fulfilled at any instance of the integration; moreover, the rela- 
tive size of the fluctuations are always kept under control during the integration meaning 
that (5f/$ <C 1 and 5v/v <C 1 are satisfied for almost all points; thus violently fluctuat- 
ing intervals where the whole procedure becomes doubtful are excluded. As was already 
anticipated in M on the classical level, we again find from the numerical integration that 
the influence of a linear potential on the dynamics during the phase of kinetic inflation 
is negligible also when fluctuations are included. The drastic deviations of the stochastic 
trajectories from the classical ones are present also for zero potential and are qualitatively 
similar. From equations (|63|) and (^) it follows that the fluctuations in v dominate over 
those in $ by one power in t as one approaches t = 0. Indeed we find from the numerics 
that the stochastic noise term in equation (^) is subdominant. 

On the basis of the intuitions associated with the standard inflationary paradigm one 
would expect the size of the fluctuations to be set by the Hubble parameter alone. How- 
ever, in our case we can see that the results depend somewhat on the parameter |$o 0) | a^ 2 . 
This is not too surprising, as there is another mass scale in the problem, namely the Planck 
scale which is actually the dynamical variable equivalent to $, hence the dependence on 
initial conditions can be seen as the dependence on the initial value of the effective Planck 
scale. Since there is no canonical way of determining a unique, internal set of initial con- 
ditions for the pre-big-bang epoch, the size of the fluctuations can only be given in terms 
of these parameters. The qualitative picture remains however unchanged no matter what 
their specific values might be, as long as they are not too exotic. From figure |3] one can 
draw the conclusion that the choice of the initial conditions merely fixes the instant of 



20 



H 




H 



-10 




Figure 3: The numerical integration of the stochastic equations of motion ( p3|) - (pq) for 
the (+) branch of the u — — 1 theory are shown (A = 1 has been chosen). In both plots the 
result of the classical equations (bold lines) are compared with the stochastic trajectories. 
The Hubble parameter is given as a function of (negative) time, (a) Various stochastic 
paths parameterized by the effective fluctuation strength |$ 



(0)1 3/2 











(as indicated) are 

compared, (b) The specific path "5.0" of (a) for different randomly chosen initializations 
of the random number generator. All stochastic trajectories emerge from the same starting 
point being also the starting point of the upper (H > 0) classical solution. 



time, when the stochastic trajectory bends away from the original classical solution. 
Finally, we should comment on the assumed value of the splitting parameter e = 1. In 
investigations of potential driven inflationary scenarios with a period very close to the 
actual de Sitter epoch one tends to go with the value of e to zero. This is usually very 
helpful in simplifying calculations, but in our opinion not so well justified on physical 
grounds, as strictly speaking in this limit all the Fourier modes except the zero mode are 
considered "fluctuations". Hence here we decided to keep e finite and of the order one 
which singles out the causal horizon as the splitting scale at any instant of the evolution. 
It is clear that e could have been chosen smaller - this would mean that we shift more and 
more degrees of freedom into the fluctuating part of the field. It turns out that over a wide 
range of finite e the results of the stochastic analysis are indeed physically equivalent. 
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5 Summary and Conclusion 



Jordan-Brans-Dicke (JBD) theories with a linearized potential for the JBD field are in- 
vestigated in the framework of a stochastic analysis which is capable of taking fluctuations 
of this field into account; their backreactions on the classical background are examined. 
We split the JBD scalar in sub- and super-horizon parts treating the long wavelength 
modes as the background quantity whose time evolution is subject to the short wavelength 
fluctuations. We derive the stochastic equations of motion for the system together with 
the two point correlation functions of the stochastic noise operators; they finally set the 
strength of the random force term in the Langevin equations which, in principle, can be 
integrated for any potential. 

For a non vanishing potential we give the scalar field mode functions for several limiting 
cases; they finally enter the expressions for the correlation functions. We compute the 
mode functions in the limit of \f\t 3> 1 for general uj and give the approximate solutions 
for y/\t <C 1 for the uj = — 1 model in terms of power law series; this allows one to 
study how the mode functions are modified if a small but finite potential is switched on. 
The mode functions can not be given analytically in the most general case; Therefore, 
we integrate the system numerically and analyze the influence of general A and to. We 
find that all mode functions are regular at t — and turn from a power law growth to an 
oscillatory regime after a time t given by the ratio |k|/A. Moreover, we demonstrate that 
they strongly depend on the value of uj chosen with a significant increase when uj drops 
even slightly below —1. 

After defining the variances of the random variables in terms of the correlation functions 
we explicitly give their form in the zero potential limit and for the stringy model uj = —1. 
Arguing that two distinct classical solutions can no longer be distinguished as soon as the 
dispersion of the fluctuations is of the same order of magnitude as the separation between 
the two classical trajectories, we show for general uj that this condition is met for all cases 
as t = is approached. We point out that this fact is not due to the mode functions (since 
they remain finite for t — *■ 0), but due to a combination of the background fields which 
gets singular close to t = 0. We find a strong enhancement of this effect for u — > 0. This is 
where many theories of gravity arising from Kaluza-Klein theories after compactification 
to four dimensions are located J|. 

Finally, for uj = —1, we perform an integration of the stochastic equations of motion 
for a finite potential demonstrating how the stochastic ensemble corresponding to the 
full quantum scalar field evolves under the influence of the scalar field fluctuations. The 
dynamics is such that the ensembles representing classical solutions which belong to dis- 
connected solutions start overlapping. We again find that the dispersion of the fluctua- 
tions grows to achieve the magnitude of the term separating the two classical solutions. 
This phenomenon can be interpreted as the quantum mechanical realization of connecting 
classically disconnected solutions at the level of field theoretical considerations. 
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